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I. INTRODUCTION 

Fifty years after Anderson's original prediction |l| of 
electron localization in weakly disordered solids, exper- 
iments demonstrated localization of ultracold neutral 
atoms in a disordered optical potential 0, Q . The great 
advantage of these atomic systems is the control of disor- 
der form and strength, interactions, dimensionality, etc., 
in principle allowing for unambiguous observation of the 
phenomenon under a variety of circumstances. However, 
the localization observed so far is qualitatively indistin- 
guishable from classical localization of particles in a dis- 
ordered trapping potential. Identifying the observed lo- 
calization as a quantum effect requires special care Q. 
It is thus of utmost interest to use the flexibility of the 
atomic systems to engineer disorder such that a unique 
signature of Anderson localization emerges. In this work 
we show both analytically and numerically that a gen- 
eral class of disorder models similar to the well-known 
random dimer model exhibits multiple derealization en- 
ergies. We propose an experimental realization of this 
class of models with ultracold atoms in optical lattices, 
using current state-of-the-art techniques. Such an ex- 
periment would feature an unusual and unmistakable ex- 
perimental signature of quantum localization in the form 
of sharp dips in the measured momentum distribution. 
The possibility of the direct observation of the quan- 
tum localization-delocalization transition in suitably de- 
signed optical lattices, particularly in situations where 
there should be no such classical transitions, is the pri- 
mary motivation of our theoretical work. 

On another front, it is very interesting to investigate 
the effect of introducing correlations to the disorder po- 
tential in Anderson- type models in higher dimensions. It 
is well-known that Anderson localization in one dimen- 
sion can be circumvented under certain conditions for 
disorder possessing long-range [4-6] or short-range 
correlations or for pseudorandom models where the dis- 
order arises from an additional potential incommensurate 



with the lattice potential [l5l-[l8|. One well-studied ex- 
ample of such one-dimensional derealization is the one- 
dimensional (ID) random dimer model [9], first studied in 
the context of the existence of extended resonance states 
in incommensurate ID lattices with complex unit cells [3] 
and later extensively applied to the metal-insulator tran- 
sition in conducting polymers [10? EH- I n this ID random 
dimer model, the localization length of the single-particle 
states diverges at a pair of critical energies. Similar be- 
havior has also been shown for a two-leg ladder version 
of the same model [l3j . In this work, we consider several 
distinct models in an N-\eg system. Such systems fill an 
interesting niche between ID and 2D derealization be- 
havior. Observing correlation-induced transitions in this 
case using the clean and controllable optical lattice sys- 
tem is an intriguing prospect. (The word "correlation" 
throughout this article and in the title of this paper refers 
exclusively to correlations in the disorder within the non- 
interacting quantum tight-binding model - mutual inter- 
action between the particles themselves, which is often 
also referred to as "correlations" in the literature, is be- 
yond the scope of the current work.) 

When one couples a set of ID chains, an additional 
degree of freedom and associated symmetries come into 
the picture. Depending on the particular experimental 
scheme, the disorder potential along the vertical (i.e., in- 
ter leg) direction may respect different symmetries. In 
the full 2D case, for instance, the potential may be invari- 
ant under the group of vertical lattice translations, T vert . 
One may then apply this potential on an N- leg lattice. 
Below we call such a potential a vertical stripe model, see 



Fig. 1(b) The vertical stripe scheme associated with the 
iV-leg model with intraleg dimer correlations turns out to 
be trivial in the sense that it can be mapped onto N sepa- 
rate dimer chains by a unitary transformation. This fact 
implies that the physical observables are essentially the 
same as those for the ID dimer model, which can also be 
inferred from the separability of the disorder potential. 

It is clear that there is an infinite number of ways in 
which this separability can be broken. To gain insight 
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FIG. 1: (Color online) Impurity pattern for a four- leg lat- 
tice with horizontal (i.e., intraleg) dimer correlations and the 
following vertical (i.e., interleg) correlations (a) uncorrelated 
model, (b) vertical stripe model, (c) diagonal stripe model, 
and (d) staggered stripe model. 



into which cases will preserve the existence of dereal- 
ization transitions and will take us out of the simple 
effectively ID dimer model, one needs to consider the 
symmetries of the model in a general case. 

In the present work we will show that breaking the 
translational invariance of the potential in the vertical 
direction and reducing the symmetry of the model down 
to T vert /Q, where Q is a discrete group, can lead to qual- 
itatively new effects that can clearly distinguish the dis- 
ordered TV-leg model from the ID random dimer model. 
We argue - and show for particular examples - that the 
physical phenomena in these schemes are determined by 
the properties of groups Q for each given number of legs, 
N. In this way one effectively characterizes properties 
of correlated disordered models based on the factoriza- 
tion of the group of translations with respect to discrete 
groups p~o| . This is one of the key results of our paper, 
and it is in some sense reminiscent of (and similar to) 
the classification of fully (uncorrelated) disordered sys- 
tems [20|, HH and topological phases [22| based on the 
Cartan classification of symmetric spaces. This is impor- 
tant especially because it is known that the type of cor- 
relations can have a strong effect on localization proper- 
ties of disordered systems. Therefore characterization of 
schemes according to their symmetry properties directly 
links these schemes with the expected universal behavior 
of observables that can be measured experimentally. Our 
current work and the earlier works on the classification 
of uncorrelated disordered systems [20| and topological 
phases [22| demonstrate that underlying group symme- 
tries provide essentially a complete classification of the 
possible universality classes of noninteracting quantum 
models although the inclusion of interaction effects in 
these group classifications remain an important unsolved 
challenge. 



In the models where the group Q = Z2 (we ter m thi s 
model a staggered vertical stripe model, see Fig. |l(d)[ ), 
the invariance under the group of translations T vert is 
broken along the vertical direction on 2 lattice spacings, 
r\n,i) = \n,i + 2), where r G T vert . This qualitatively 
changes the physical properties of the system compared 
to the vertical stripe model, where the translational in- 
variance in the vertical direction holds for all translations. 
Breaking of the translational invariance (i) reduces the 
symmetry of the model from T vert down to T vert /Z2] (ii) 
preserves the existence of truly delocalized states in a 
model with M-mers, as we will see from our numerical 
calculations of the localization length for the staggered 
vertical stripe model; and (hi) leads to qualitatively dif- 
ferent forms of the momentum distributions for different 
number of legs, TV, and correlations, M. 

One could also generalize the staggered model and con- 
sider, for example, longer periodic correlations of the site 
energies, e n ^, in the vertical direction. For example one 
can violate the translational invariance in the vertical di- 
rection on three and more (m in general) lattice spacings, 
which will reduce the symmetry by the discrete group 
Z m down to the factor space T vert /Z m . As in the ver- 
tical stripe model case, the presence of this symmetry 
is the origin of the appearance of delocalized states in a 
corresponding M-mer (M = 2, 3, ..) model. 

In this work we will discuss in detail and contrast two 
simple but spectacular cases out of the infinite number of 
ways in which the separability of the disorder potential 
can be broken. In the first example, w hich is termed the 
N-leg diagonal stripe model (see Fig. |l(c) ), the separa- 
bility is broken but also the invariance under the group 
of translations T vert in the vertical direction is broken 
completely (as opposed to the case where it is broken 
down to Tvert/Zm). We show that such a modification 
leads to the complete localization of states. The sec- 
ond example is the TV-leg staggered vertical stripe model 
which clearly yields new experimentally observable ef- 
fects which we discuss in detail. 

Generally, violation of the translational invariance in 
the vertical direction by any discrete group, Q, will work 
in a similar way and one should expect the presence of 
delocalized states here as well. Amazingly, this property 
implies the classification of the correlated disordered sys- 
tems characterized by the orientation-preserving trans- 
formations corresponding to Q. 

The paper is organized as follows. First, we pro- 
pose a procedure to experimentally realize the ID ran- 
dom dimer model in an optical lattice and measure the 
derealization energies in Sec. [Ill This model offers a 
unique opportunity to both explore correlation-induced 
derealization and observe an unambiguous signal of An- 
derson localization. Second, in Sec. IIIII we generalize 
to anisotropic short-range-correlated disorder in an TV- 
leg ladder. We analytically and numerically investigate 
localization properties and propose experimental tech- 
niques to realize these models with ultracold atoms in 
an optical lattice. We conclude in Sec. [TV] with a brief 
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summary. Exact solution of the vertical stripe model and 
some technical details are given in the Appendices. 



II. ID RANDOM DIMER MODEL 

For a ID quantum system with uncorrelated Anderson 
disorder, the wavefunction is always localized at any en- 
ergy and behaves asymptotically as ipE (z) ~ e~ z ^^ E \ 
where the localization length, £ generally depends 
on the energy. In cold atom systems, the density pro- 
file is a natural quantity to measure, so the localization 
length can be measured as a function of disorder 0, Q , 
although in practice all measurements are averages over 
energy due to the distribution of the initial wavepacket 
[23|,|23. However, a classical particle in a disordered po- 
tential can also be localized simply by random potential 
peaks higher than the particle energy. The localization 
length of a classical particle increases with energy, as a 
larger region becomes classically accessible, but the same 
is true for a quantum particle. Therefore, without suffi- 
cient knowledge of the potential landscape and the par- 
ticle energy, it is difficult to determine whether or not 
the localization seen in a given cold atom experiment is 
a quantum effect [25j. 

However, there are some correlated forms of disorder 
for which the ID Anderson localization length does not 
depend monotonically on energy. The most famous of 
these is the random dimer model [|[, 

L 

H = Yl [ 6 ^ C n C n + t {c ] n C n+1 + 4 +1 C n ) , (1) 
n=l 

where L is the length of the chain, t is the hopping pa- 
rameter, and e n = ±e is the random on-site potential. In 
this model the binary disorder has short-range correla- 
tions such that equal energy sites always occur in pairs. 

Although uncorrelated disorder localizes all ID states 
to some finite length, in the random dimer model the 
dimer correlation allows a pair of transmission resonances 
at critical energies E c = d=e 04H1- There is a boundary 
value, e = t, such that for larger e there are no delocalized 
states [lj, LL3| , but for e < t (the "bulk" region) , the lo- 
calization length actually diverges as £ (E) ~ ( E ~ t Ec ) 
Right at the boundary, e = t, 

Thus the localization length has peaks, as shown in 
Fig. [2j and we see that for correlated disorder the local- 
ization length can actually decrease with energy in some 
cases. Clearly this is a purely quantum signature as- 
sociated with true Anderson localization as opposed to 
classical localization, as noted in Ref. [25j where certain 
long-range disorder correlations are predicted to give rise 
to a qualitatively similar phenomenon. However, in the 
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FIG. 2: (Color online) Localization length versus energy for 
random dimer disorder. Here e = t/2, a is the lattice spacing, 
and the two values of the on-site energy occur with equal 
probability. 



case of Ref. [25| it is not easy to directly measure local- 
ization length as a function of energy because the atoms 
have a broad initial distribution of energies [24]. For 
the random dimer model, though, the existence of de- 
localization points allows the anomalous behavior of the 
correlation length to be inferred easily: After releasing 
the atoms from the initial confinement, one simply waits 
long enough that the atoms with energies near the crit- 
ical derealization energies have propagated out of the 
lattice, so that the final momentum distribution for the 
remaining atoms has dips which will be evident in a time- 
of- flight image. 

We now turn to the details of realization and ob- 
servation of the above physics in optical lattice experi- 
ments. Disordered potentials can be realized in a cold 
atom system either by directly disordering the optical 
trap II or by loading a regular lattice with one mo- 
bile atomic species and one immobile atomic species that 
serves as a pinned impurity distribution for the mobile 
atoms [3, HH, [27} . The controllability of nearly all as- 
pects of optically trapped cold atomic systems allows sev- 
eral complementary schemes to realize disordered poten- 
tials. Below we outline one possibility, sketched in Fig. [3] 

To implement ID random dimer model disorder, one 
could pin impurity atoms in a ID optical lattice. Stan- 
dard techniques |28i-l3Qj allow selective removal of atoms 
such that only empty or doubly-occupied sites remain, 
with a random distribution. Each well is then split by 
adiabatically turning off the original lattice while turning 
on a lattice with half its wavelength [3l| . Randomly dis- 
tributed pairs of singly occupied sites and pairs of empty 
sites result if the impurity atoms are repulsively inter- 
acting or if a second round of purification is performed 
to remove doubly occupied sites. (A similar procedure 
has recently been proposed to realize the dual random 
dimer model [14], where sites with +e never occur in 
pairs, but that proposal introduces unwanted additional 
correlations.) Errors in the disorder initialization may re- 
sult in additional uncorrelated disorder with a mean free 
path i. The effect of these errors will be to superimpose 




FIG. 3: (Color online) Scheme to realize random dimer disor- 
der, (a) Quench to random static distribution, (b) Purify to 
empty or doubly occupied sites, (c) Split wells adiabatically. 
(d) Random dimer disorder. 

another overall localization envelope on the wavefunc- 
tion, e~ z ^ . We will see that this is not a handicap, as 
long as £ is not small compared to the system size. 

Once the disordered potential is produced, the mobile 
atoms - bosons with intraspecies scattering length tuned 
to zero - can be loaded into the lattice in the presence of 
an additional tightly confining potential, as in Refs. [HQ- 
The final distribution of crystal momenta for the local- 
ized atoms to be imaged is 

r L/2 

p f (K z ) = Pl (K z ) / dz\^ Kz {z)\ 2 

J-L/2 

= Pi (K z ) J dE5(E -2t cos K z a) 
x [ L/2 dz (1/e (E) + l/£) e-^/^+i/') 

J-L/2 

= (l - e -L(imtcosK z a) + l/i)^ p . {Rz) (3) 

Here pi (K z ) is the initial momentum distribution, and 
we have assumed only the first band is populated. Adia- 
batically turning off the lattice maps quasi- momentum to 
real momentum (3(J [32| and taking a time-of- flight image 
of these atoms (after turning off interaction with impu- 
rity atoms, if applicable) yields a momentum distribution 
with marked dips at momenta 

k z = a -1 arccos (±e/2£) , (4) 

corresponding to the derealization energies. We show an 
example of this in Fig. |4] for an initial gaussian momen- 
tum distribution. Note that this signature is somewhat 
suppressed by any non-dimer impurities present, but is 
still visible as long as the system size is not much larger 
than the mean free path of the non-dimer impurities. The 
dips also become narrower for larger systems. These dips 
are an unmistakeable signature of quantum localization- 
delocalization points. 
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FIG. 4: (Color online) Momentum distribution along the z- 
axis of atoms in random dimer disorder with e = t/2. The 
initial momentum distribution is taken here to be gaussian. 
The distribution measured via time-of-flight after the delocal- 
ized atoms have left the system is shown in the presence of 
additional non-dimer impurities with mean free path i. The 
system size is taken to be 100 sites. 



III. N-LEG ANISOTROPIC ALLY-CORRELATED 
DISORDER 

We now generalize to an Anderson disordered system 
on an A 7 "- leg optical lattice, which may be thought of as a 
type of intermediate system between ID and 2D. We will 
again examine the role of correlations in the disorder po- 
tential. We consider an extreme case of anisotropic cor- 
relations with short-range dimer correlations along each 
leg and infinite-range correlations between legs. 

The Hamiltonian for this system is 

L N , 

71=1 2=1 ^ 

+ t±(4i,i C n,i+l + C l,i+l C n,i)|. (5) 

Here t and t± are constant intraleg and interleg hop- 
ping parameters, respectively, L is the number of sites 
in each leg, and N is the number of legs in the ladder. 
Energy e n ^ = ±e is the random on-site potential. De- 
noting by ?p nj i the single-particle wave function at site n 
of the i-th leg and defining the two-column wavefunction 
C0n,i, • • • ,i/>n,N, ^n-i,i, • • • , ^n-i,iv) , one can write 
the solution of the Schroedinger equation H\ip) = E\i/j), 
where |^) = En^VwK^ as ^™+i = u where 

T{ is the transfer matrix. 

Correlations in the random on-site potential are im- 
posed along the length of the chain such that equal en- 
ergy sites occur in pairs, or more generally in "M-mers" . 
Correlations are imposed between chains in various ways 
to be discussed below. 
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A. Vertical stripe model 

The simplest example of cross-leg correlations is a sys- 
tem where e n ^ = e n , i.e., all chains are identical copies of 
a ID random dimer potential. In this case, the transfer 
matrix is 



T 

-J- n 



t t 



1 + CT3 



iI N ®a 2 , (6) 



where cr 12 ,3 are the Pauli matrices, In is an N x N 
identity matrix, and S is an TV- dimensional matrix re- 
sponsible for transport in the vertical direction with 

One can realize this system as in the ID case, except 
that now one must initialize randomly distributed stripes 
of empty sites and doubly occupied sites. This could be 
done by loading the lattice with two atoms on every site 
[29j and temporarily turning on a optical speckle poten- 
tial that is disordered along the length of the chains but 
constant across the legs, as in Ref. |2[ (see Fig. 5(a) ). The 
speckle beam introduces a differential light shift so that 
each cross-leg column of sites has an associated random 
atomic transition frequency. By slowly sweeping the fre- 
quency of an external microwave field over an appropri- 
ate range, all atoms whose internal transition frequency 
is shifted into a given range by the speckle will be adi- 
abatically transferre d into a different internal state, as 
in Ref. [33j (see Fig. |5(b)[ ). Upon adiabatically turning 
off the microwave field and the speckle, one is left with 
randomly distributed stripes of sites with atoms in a dif- 
ferent internal state. These atoms can be removed with 
a blast pulse. Splitting each well in the intrachain direc- 
tion as before results in the desired disorder. We sketch 
the process in Fig. [5] 

Note that, for a continuous distribution of random 
speckle shifts, there will always be a fraction of columns 
for which the adiabatic sweep results in a mixed state, 
as shown in Fig. 5(b) This leads to imperfections in the 
disorder pattern after the projective blast pulse. If one 
sweeps the microwave field 60 kHz over 20 ms with pulse 
parameters as in Ref. (H], the detuning range that gets 
transferred to a mixed state is roughly 2 kHz wide. As- 
suming the random light shifts have a frequency range 
similar to that of the detuning sweep (corresponding to 
a 2 mW speckle beam [34| in this case), and accounting 
for errors from nonadiabatic transitions as well, one ob- 
tains an average distance between errors in the random 
dimer disorder of i ~ 70 lattice spacings. However, t can 
be greatly increased by using a slower sweep. 

Since the potential is separable in the inter- and intra- 
chain directions, one simply obtains the same localization 
length as a function of energy along the chain as in the 
ID case above (see App. [A]), and the same dips in the 
momentum profile. The only difference is that there will 
be N times as many delocalized states, due to the N 
possible transverse modes. 
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FIG. 5: (Color online) Scheme to realize iV-leg random dimer 
disorder, (a) Load n — 2 Mott insulator and apply speckle 
beam, (b) Sweep detuning to adiabatically transfer atoms 
in randomly selected columns from internal state \a) to \b). 
State transfer occurs for atoms in sites that are shifted into 
a window of detunings near resonance by the speckle beam. 
The internal state is mixed for a small window of intermediate 
detunings. (c) Remove atoms in state \b) [blue (dark gray)], 
(d) Split wells adiabatically, resulting in dimer disorder. 



B. Diagonal stripe model 

Clearly, we must break the separability of the disorder 
potential in the intra- and interchain directions to obtain 
different behavior. The simplest way to do that is to 
rotate the propagation axis of the speckle beam in the 
2D plane such that it is no longer perpendicular to the 
chains. For instance, if the same procedure as above is 
performed with the speckle beam pointed at a 45° angle 
to the legs in the plane, then each leg will have identical 
dimer sequences (except for at the boundaries), but the 
sequence on each leg will be shifted by two lattice sites 
relative to the leg above it, such that e n ^ = e n _2,i+i, i = 
1, • • • N — 1. Here, the transfer matrix is given by 



N 



t ' 



1 + CT3 



N 



><T2, (7) 



where Y n ^ is an A/"-dimensional diagonal matrix with 



entries: (Y n ^)j 



i+2(j-l) 



We have used this form of the transfer matrix to cal- 
culate the localization length of a two-leg dimer system 
numerically; see Fig. [6j In Fig. [6fa) we see the presence 
of derealization transition in this model at zero-energy 
when e = t±. This point is special in the sense that it 
corresponds to the case when the on site chemical poten- 
tial is sufficiently large to guarantee for a particle to hop 
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FIG. 6: (Color online) Localization lengths for two-leg dimer 
(N = 2, M — 2) diagonal stripe model at: (a) e = t/2 and 
t ± t/2; (b) e = t/2 and £_l = t; (c) e = 3t/4 and t± = t. 
Numerical calculation was done for a system with 8000 sites 
in the horizontal direction. 



C. Staggered vertical stripe model 

Another interesting disorder pattern that is not sep- 
arable is dimer-correlated disorder in each chain with 
long-range interchain correlations such that each chain 
is the mirror image of the one above it, i.e., the disor- 
der potential is staggered along the interchain links. The 
Hamiltonian of the model is defined by Eq. (J5j) where the 
site energies e n ^ = (— l) l e n along the vertical direction 
i = 1, ..TV depend on a single random parameter e n = ±e. 
The transfer matrix of the model reads 



EI 



N 



eJ N 



t 



t 



1 + CT3 



UN'S) cr 2 , (8) 



where In is identity matrix, while (JN)jj/ — (— 

In Fig. we present data for localization lengths in a 
staggered model for two- leg dimer (a) and trimer (b), and 
three-leg dimer (c) and trimer (d) models, respectively. 
All numerical simulations here have been done for a chain 
of 6000 M-mers. We see here a very rich collection of 
true delocalized states. It is important to notice that the 
staggered model can not be reduced to a set of N one-leg 
models as in the plain vertical stripe case. 

Fig. [7] clearly shows the presence of derealization of 
the states in all cases (a-d), however another interest- 
ing property of the staggered model is evident here. 
Namely, we see a region of well-extended states around 
E ~ 1.5t — 3t. Their localization lengths are comparable 
with the length of the chain under consideration. This 
fact produces large regions of dips in a momentum dis- 
tribution of the states (see Fig. [8j). 

Fig. [8] represents distributions of horizontal (i.e., in- 
terleg) momentum after delocalized and extended states 
have left a sample with 100 sites. We take an initially 
Gaussian momentum distribution centered at zero. 

From Appendix [Bj we have that the connection of the 
derealization energy E c and the derealization momen- 
tum p c is defined by the expression 



arccos 



e c ± 

2t 



(9) 



for N = 2, and 



freely in the vertical direction. This effectively restores 
the translational invariance in the vertical direction and 
makes the model analogous to the vertical stripe model. 
Supporting evidence is that the critical index correspond- 
ing to this transition is v = 2, which coincides with that 
of ID dimer model in the "bulk" region. For other val- 
ues of model parameters the single particle states remain 
localized; see Fig. [6f b) and (c) for two particular choices 
of model parameters. 



Pc 



L 2t J 



(10) 



2t 



for N = 3. It is important to note here that the spectrum 
depends on the number of legs, A 7 ", in a nontrivial way. 
This implies Af-dependence of observables in general and 
of the momentum distribution in particular (see Fig. [8]). 

A staggered disorder potential can be realized experi- 
mentally by loading the lattice into a n = 1 Mott insula- 
tor state, then using the speckle beam to select columns 
of sites to flip to a different internal state, as in the verti- 
cal stripe case discussed above. However, here one must 
specify the two relevant internal states to have opposite 



7 



f(E)/fl 




P(*) 



«E)/fl 



o 

E/t 



12000 




(c) 




f(E)/a 



-3 







3 






(d) 


K 




M 


03D - 


-3 







3 



FIG. 7: (Color online) Localization lengths of staggered ver- 
tical stripe models for e = t/2. (a) Two- leg dimer model 
(TV = 2, M = 2) showing E c ~ 0, (b) two- leg trimer model 
(TV = 2, M = 3) showing E c ~ ±0.12t, (c) three-leg dimer 
model (JV = 3, M = 2) showing £ c ~ ±0.5*. (d) two-leg 
trimer model (N = 2, M = 3) showing £ c ~ ±0.5£. 
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FIG. 8: (Color online) Momentum distribution after delocal- 
ized states have left the sample for N=2, M=2 (a); N=2, 
M=3 (b); and N=3, M=2 (c) cases, respectively. The initial 
distribution is the same as in Fig. [H 



Zeeman shifts. Then, deforming the lattice to a two- 
period superlattice, ramping up the tunneling, applying a 
magnetic field gradient, and ramping the tunneling back 
down results in the atoms being shifted according to their 
internal state, as shown in Fig. [9j This is somewhat sim- 
ilar to the procedure used to initialize a Neel state in 
Ref. [il. 

More straightforward methods of realizing the short- 
range correlated disorder are also available, though they 
require sophisticated single-site control. In that case, 
one can create completely arbitrary disorder patterns. 
One pos sibility is to use a custom holographic mask as in 
Ref. [36] to generate a two-dimensional (2D) lattice po- 
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FIG. 9: (Color online) Scheme to realize staggered iV-leg 
dimer disorder, (a) Load n — 1 Mott insulator, (b) Flip inter- 
nal state of atoms in randomly selected columns [blue (dark 
gray)] as in Fig. [5] and lower every other interleg barrier, (c) 
Apply magnetic field gradient to obtain state-dependent tilt, 
(d) Restore barriers, flip back internal states, and split wells 
adiabatically. 



tential with the desired disorder pattern. Another is to 
use a regular 2D optical lattice, but to etch the desired 
disorder pattern onto a regular array of impurity atoms 
with single-site resolution [33|, l37j . 



IV. CONCLUSIONS 

We discussed a class of quasi- ID disordered models 
with correlations between random variables, which, de- 
pending on the symmetries of the model which we ex- 
plicitly identify and discuss, may or may not exhibit 
localization-delocalization transitions. We establish us- 
ing a few examples how such models and their observ- 
able characteristics (such as momentum distribution and 
single particle spectrum) strongly depend on symmetry 
properties of the effective hopping Hamiltonian, which, 
in turn, classify the models accordingly. 

We have shown that cold atoms trapped in an optical 
lattice can realize the well-known one-dimensional ran- 
dom dimer model, which exhibits critical derealization 
energies, and proposed an experimental procedure to ob- 
serve the derealization via time-of-flight imaging. We 
have also generalized the random dimer model to TV- 
leg systems with highly anisotropic correlated disorder 
and found new transmission resonances induced by the 
transverse degree of freedom. We have calculated the 
localization length as a function of energy and identi- 
fied the critical energies. These resonances, both in the 
ID and TV- leg systems, present unique and unambiguous 
signatures of the quantum derealization energies. Fur- 
thermore, the iV-leg case provides an exciting prospect 



for experiments to explore correlation-induced derealiza- 
tion outside of the strictly ID limit with the current state 
of the art. Also, any experimental observation of criti- 
cal energies separating localized and delocalized states in 
well-designed disordered optical lattices along the line of 
our suggestions in this work will immediately establish 
the occurrence of the quantum localization phenomena 
in cold atomic sysems. 

We thank T. Porto and S. Rolston for helpful dis- 
cussion. This work supported by AFOSR-MURI, ARO- 
MURI, DARPA-OLE, and NSF-JQI-PFC. 



Appendix A 

Localization properties of the system can be investi- 
gated by calculating Lyapunov exponents. According to 
Oceledec's theorem [38|, the power of the product of 
L random transfer matrices has a limiting value 



A = lim 



-i 1/2L 



n=l 



(Al) 



The closest to unity eigenvalue of A can be represented 
as e x ^ E \ where \{E) is the Lyapunov exponent and the 
localization length is simply its inverse, £(E) = 1/X(E). 
A derealization transition occurs at energies E c where 
the localization length diverges. We analytically calcu- 
late the critical energies, E C1 for iV-leg M-mer disorder 
in two steps. First we reduce the iV-leg problem to an 
effective one-leg problem. Then we solve the problem for 
arbitrary M-mers in a one-leg chain. 
The eigenvalues of t±S are 



2t_L cos 



E 



( irk 



N,k 



t± cos (irk) , 




k = 1,2 



, N, N>2 

N = 2 
N = 1 



so the whole transfer matrix can be recast in the form 
of N independent ID transfer matrices corresponding to 



(k) 

different transverse modes, T n = @ fe Tn , where 



T (k) 



<E 



EN,k 



t 



a 3 



KJ2- 



(A2) 



Therefore, the critical energies, E c , of the ladder model 
correspond to those of a single decoupled chain plus the 
transverse mode energy. 

The M-mers are now treated by defining a new basic 
transfer matrix, which is the Mth power of the expression 
(|A2p , [Tn^] , with randomly distributed site energies 
e n . At a derealization energy, E C1 the transfer matrix 
has an eigenvalue being ±1 [7] (in units of t). Here this 
implies that for one of the values of fc, 



[T«] M = ±J 2 



(A3) 
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To solve this equation it is convenient to express the 
transfer matrix ()A2|) in terms of trigonometric functions, 

T (k) = e i«.o- = h CQSUJ + s[nuJ ^ (A4) 

where oj = arccos ^ E - E ^k-e n ^ an( ^ ^ _ 

( '-ii^J'-^ cot ^)- 

Then [T^ k) ] M = I 2 cos (Mw) + z^f sin (Mw). There- 
fore, the delocalization condition is satisfied for cj = 
7rg/M, g = 1, . . . , 2M, and the critical energies are 

7T<7 

£ c 2t cos — + £w >fe ± e, 
g = 1, . . . , 2M - 1; q ^ M, and fc = 1, . . . , TV. (A5) 

Note that the trivial solutions uo = 0,7r,27r, correspond- 
ing to transfer matrix T n = ±l"2iv, should not be con- 
sidered. Below we write explicitly the critical points for 
dimers, trimers, and tetramers: 

M = 2 : E c = E Njk ± e; 

M = 3 : £ c = ± £_l ± e; (A6) 
M = 4: E c = ± >/2*_l ± e, £;v,fc±e. 

Generally speaking, Eq. (|A3|) is a sufficient condition 
for having a delocalization transition but it is not a nec- 
essary one. Therefore, we have checked numerically for 
up to TV = 7 legs and for correlations of up to M = 5 
sites that Eq. (|A5[) captures all the critical points. 

Away from the critical points, the localization length 
must be calculated numerically from Eq. (|A1[) . Fig. [2] 
shows some numerical results for localization lengths in 
a chain of L = 5000 sites. The critical energies obtained 
are in perfect agreement with the analytical result of 
Eq. ((M]) - 



model defined by the Hamiltonian (j5j) with e n ^ = (— l)*e 
(with non-random e). The reason for this is the follow- 
ing. Theoretically we know that the particles in a disor- 
dered chain are either localized (£(E) is finite) or delo- 
calized (£(E C ) —> oo). Experimentally the localized on- 
shell particles are being "detected" by measuring their 
momentum distribution profile. At the time when the 
measurement is done the energy of particles is linked to 
the momentum in a way which has no information about 
disorder and its pre-history, and therefore the spectrum 
of the detected particle precisely coincides with the spec- 
trum of the free particle in a non-random model. 

As the first step we Fourier transform the Hamiltonian 
with no disorder in the horizontal direction and intro- 
duce dimensionless momenta, p z . Then the Hamiltonian 
for the two leg dimer system, N = 2, (eigenfunctions of 
which are ipi iP , i = 1, 2) becomes 

H= (t ± -e) + 2 * cos b*]l- (Bl) 

Diagonalization of this Hamiltonian yields the following 
dispersion relation 

E = ±^/e 2 -Ki + 2£cos[^] (B2) 

which produces Eq. (j9j). 

For three leg system we have the following Hamiltonian 
in the momentum space 

/ c t± \ 
H= t± -e £_l +2tcos[p 2 ]i. (B3) 



Appendix B 

In this Appendix we will calculate the spectrum of a 
particle in the staggered vertical stripe model (class Z2). 

In order to obtain the relation between the energy and 
the momentum of the free on-shell particle in a stag- 
gered disordered model, one needs to calculate the disper- 
sion of the free particle in the staggered but non-random 



Upon diagonalization of Eq. (jB3|) we arrive to the follow- 
ing spectrum 

E _ f ± ^/2ii + <? + 2t cos[p z ] , B4 . 
which defines Eq. ([TO]) of the main text. 
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